home *** CD-ROM | disk | FTP | other *** search
/ Aminet 21 / Aminet 21 (1997)(GTI - Schatztruhe)[!][Oct 1997].iso / Aminet / gfx / show / gs_src_gs.lha / gs5.03 / spdiff.c < prev    next >
C/C++ Source or Header  |  1996-04-18  |  8KB  |  298 lines
  1. /* Copyright (C) 1994, 1995, 1996 Aladdin Enterprises.  All rights reserved.
  2.   
  3.   This file is part of Aladdin Ghostscript.
  4.   
  5.   Aladdin Ghostscript is distributed with NO WARRANTY OF ANY KIND.  No author
  6.   or distributor accepts any responsibility for the consequences of using it,
  7.   or for whether it serves any particular purpose or works at all, unless he
  8.   or she says so in writing.  Refer to the Aladdin Ghostscript Free Public
  9.   License (the "License") for full details.
  10.   
  11.   Every copy of Aladdin Ghostscript must include a copy of the License,
  12.   normally in a plain ASCII text file named PUBLIC.  The License grants you
  13.   the right to copy, modify and redistribute Aladdin Ghostscript, but only
  14.   under certain conditions described in the License.  Among other things, the
  15.   License requires that the copyright notice and this notice be preserved on
  16.   all copies.
  17. */
  18.  
  19. /* spdiff.c */
  20. /* Pixel differencing filters */
  21. #include "stdio_.h"        /* should be std.h, but needs NULL */
  22. #include "strimpl.h"
  23. #include "spdiffx.h"
  24.  
  25. /* ------ PixelDifferenceEncode/Decode ------ */
  26.  
  27. private_st_PDiff_state();
  28.  
  29. #define ss ((stream_PDiff_state *)st)
  30.  
  31. /* Define values for case dispatch. */
  32. #define cBits1 0
  33. #define cBits2 4
  34. #define cBits4 8
  35. #define cBits8 12
  36. #define cEncode -1
  37. #define cDecode 15
  38.  
  39. /* Set defaults */
  40. private void
  41. s_PDiff_set_defaults(stream_state *st)
  42. {    s_PDiff_set_defaults_inline(ss);
  43. }
  44.  
  45. /* Common (re)initialization. */
  46. private int
  47. s_PDiff_reinit(stream_state *st)
  48. {    ss->row_left = 0;
  49.     return 0;
  50. }
  51.  
  52. /* Initialize PixelDifferenceEncode filter. */
  53. private int
  54. s_PDiffE_init(stream_state *st)
  55. {    long bits_per_row =
  56.         ss->Colors * ss->BitsPerComponent * (long)ss->Columns;
  57.     static const byte cb_values[] =
  58.         { 0, cBits1, cBits2, 0, cBits4, 0, 0, 0, cBits8 };
  59.     ss->row_count = (uint)((bits_per_row + 7) >> 3);
  60.     ss->end_mask = (1 << (-bits_per_row & 7)) - 1;
  61.     ss->case_index =
  62.         cb_values[ss->BitsPerComponent] + ss->Colors + cEncode;
  63.     return s_PDiff_reinit(st);
  64. }
  65.  
  66. /* Initialize PixelDifferenceDecode filter. */
  67. private int
  68. s_PDiffD_init(stream_state *st)
  69. {    s_PDiffE_init(st);
  70.     ss->case_index += cDecode - cEncode;
  71.     return 0;
  72. }
  73.  
  74. /* Process a buffer.  Note that this handles both Encode and Decode. */
  75. private int
  76. s_PDiff_process(stream_state *st, stream_cursor_read *pr,
  77.   stream_cursor_write *pw, bool last)
  78. {    register const byte *p = pr->ptr;
  79.     register byte *q = pw->ptr;
  80.     int rcount, wcount;
  81.     register int count;
  82.     int status = 0;
  83.     register byte s0 = ss->s0;
  84.     register byte t;
  85.     byte save_last;
  86.     const byte end_mask = ss->end_mask;
  87.  
  88. row:    if ( ss->row_left == 0 )
  89.       ss->row_left = ss->row_count,
  90.       s0 = ss->s1 = ss->s2 = ss->s3 = 0;
  91.     rcount = pr->limit - p;
  92.     wcount = pw->limit - q;
  93.     if ( ss->row_left < rcount )
  94.       rcount = ss->row_left;
  95.     count = (wcount < rcount ? (status = 1, wcount) : rcount);
  96.     ss->row_left -= count;
  97.     if ( ss->row_left == count )
  98.       save_last = p[count];
  99.  
  100.     /*
  101.      * Encoding and decoding are fundamentally different.
  102.      * Encoding computes E[i] = D[i] - D[i-1];
  103.      * decoding computes D[i] = E[i] + D[i-1].
  104.      * Nevertheless, the loop structures are similar enough that
  105.      * we put the code for both functions in the same place.
  106.      */
  107.  
  108. #define loopn(n, body)\
  109.   while ( count >= n ) p += n, q += n, body, count -= n
  110.  
  111.     switch ( ss->case_index )
  112.     {
  113.  
  114.     /* 1 bit per component */
  115.  
  116. #define eloop1(ee)\
  117.   loopn(1, (t = *p, *q = ee, s0 = t))
  118.  
  119.     case cEncode+cBits1+1:
  120.         eloop1(t ^ ((s0 << 7) | (t >> 1)));
  121.     case cEncode+cBits1+2:
  122.         eloop1(t ^ ((s0 << 6) | (t >> 2)));
  123.     case cEncode+cBits1+3:
  124.         eloop1(t ^ ((s0 << 5) | (t >> 3)));
  125.     case cEncode+cBits1+4:
  126.         eloop1(t ^ ((s0 << 4) | (t >> 4)));
  127.  
  128. #define dloop1(te, de)\
  129.   loopn(1, (t = te, s0 = *q = de)); break
  130.  
  131.     case cDecode+cBits1+1:
  132.         dloop1(*p ^ (s0 << 7),
  133.             (t ^= t >> 1, t ^= t >> 2, t ^ (t >> 4)));
  134.     case cDecode+cBits1+2:
  135.         dloop1(*p ^ (s0 << 6),
  136.             (t ^= (t >> 2), t ^ (t >> 4)));
  137.     case cDecode+cBits1+3:
  138.         dloop1(*p ^ (s0 << 5),
  139.             t ^ (t >> 3) ^ (t >> 6));
  140.     case cDecode+cBits1+4:
  141.         dloop1(*p ^ (s0 << 4),
  142.             t ^ (t >> 4));
  143.  
  144.     /* 2 bits per component */
  145.  
  146. #define add4x2(a, b) ( (((a) & (b) & 0x55) << 1) ^ (a) ^ (b) )
  147. #define sub4x2(a, b) ( ((~(a) & (b) & 0x55) << 1) ^ (a) ^ (b) )
  148.  
  149.     case cEncode+cBits2+1:
  150.         eloop1((s0 = (s0 << 6) | (t >> 2), sub4x2(t, s0)));
  151.     case cEncode+cBits2+2:
  152.         eloop1((s0 = (s0 << 4) | (t >> 4), sub4x2(t, s0)));
  153.     case cEncode+cBits2+3:
  154.         eloop1((s0 = (s0 << 2) | (t >> 6), sub4x2(t, s0)));
  155.     case cEncode+cBits2+4:
  156.         eloop1(sub4x2(t, s0));
  157.  
  158.     case cDecode+cBits2+1:
  159.         dloop1(*p + (s0 << 6),
  160.             (t = add4x2(t >> 2, t),
  161.             add4x2(t >> 4, t)));
  162.     case cDecode+cBits2+2:
  163.         dloop1(*p, (t = add4x2(t, s0 << 4),
  164.                add4x2(t >> 4, t)));
  165.     case cDecode+cBits2+3:
  166.         dloop1(*p, (t = add4x2(t, s0 << 2),
  167.                add4x2(t >> 6, t)));
  168.     case cDecode+cBits2+4:
  169.         dloop1(*p, add4x2(t, s0));
  170.  
  171. #undef add4x2
  172. #undef sub4x2
  173.  
  174.     /* 4 bits per component */
  175.  
  176. #define add2x4(a, b) ( (((a) + (b)) & 0xf) + ((a) & 0xf0) + ((b) & 0xf0) )
  177. #define add2x4r4(a) ( (((a) + ((a) >> 4)) & 0xf) + ((a) & 0xf0) )
  178. #define sub2x4(a, b) ( (((a) - (b)) & 0xf) + ((a) & 0xf0) - ((b) & 0xf0) )
  179. #define sub2x4r4(a) ( (((a) - ((a) >> 4)) & 0xf) + ((a) & 0xf0) )
  180.  
  181.     case cEncode+cBits4+1:
  182.         eloop1( ((t - (s0 << 4)) & 0xf0) | ((t - (t >> 4)) & 0xf) );
  183.     case cEncode+cBits4+2:
  184.         eloop1(sub2x4(t, s0));
  185.     case cEncode+cBits4+3:
  186.     {    register byte s1 = ss->s1;
  187.         loopn(1, (t = *p,
  188.               *q =
  189.                 ((t - (s0 << 4)) & 0xf0) | ((t - (s1 >> 4)) & 0xf),
  190.               s0 = s1, s1 = t));
  191.         ss->s1 = s1;
  192.     }    break;
  193.     case cEncode+cBits4+4:
  194.     {    register byte s1 = ss->s1;
  195.         loopn(2,
  196.               (t = p[-1], q[-1] = sub2x4(t, s0), s0 = t,
  197.                t = *p, *q = sub2x4(t, s1), s1 = t));
  198.         ss->s1 = s1;
  199.     }    break;
  200.  
  201.     case cDecode+cBits4+1:
  202.         dloop1(*p + (s0 << 4), add2x4r4(t));
  203.     case cDecode+cBits4+2:
  204.         dloop1(*p, add2x4(t, s0));
  205.     case cDecode+cBits4+3:
  206.     {    register byte s1 = ss->s1;
  207.         loopn(1, (t = (s0 << 4) + (s1 >> 4),
  208.               s0 = s1, s1 = *q = add2x4(*p, t)));
  209.         ss->s1 = s1;
  210.     }    break;
  211.     case cDecode+cBits4+4:
  212.     {    register byte s1 = ss->s1;
  213.         loopn(2,
  214.               (t = p[-1], s0 = q[-1] = add2x4(s0, t),
  215.                t = *p, s1 = *q = add2x4(s1, t)));
  216.         ss->s1 = s1;
  217.     }    break;
  218.  
  219. #undef add2x4
  220. #undef add2x4r4
  221. #undef sub2x4
  222. #undef sub2x4r4
  223.  
  224.     /* 8 bits per component */
  225.  
  226. #define encode8(s, d) (q[d] = p[d] - s, s = p[d])
  227. #define decode8(s, d) s = q[d] = s + p[d]
  228.  
  229.     case cEncode+cBits8+1:
  230.         loopn(1, encode8(s0, 0));
  231.         break;
  232.     case cDecode+cBits8+1:
  233.         loopn(1, decode8(s0, 0));
  234.         break;
  235.     case cEncode+cBits8+2:
  236.     {    register byte s1 = ss->s1;
  237.         loopn(2, (encode8(s0, -1), encode8(s1, 0)));
  238.         ss->s1 = s1;
  239.     }    break;
  240.     case cDecode+cBits8+2:
  241.     {    register byte s1 = ss->s1;
  242.         loopn(2, (decode8(s0, -1), decode8(s1, 0)));
  243.         ss->s1 = s1;
  244.     }    break;
  245.     case cEncode+cBits8+3:
  246.     {    register byte s1 = ss->s1, s2 = ss->s2;
  247.         loopn(3, (encode8(s0, -2), encode8(s1, -1),
  248.               encode8(s2, 0)));
  249.         ss->s1 = s1, ss->s2 = s2;
  250.     }    break;
  251.     case cDecode+cBits8+3:
  252.     {    register byte s1 = ss->s1, s2 = ss->s2;
  253.         loopn(3, (decode8(s0, -2), decode8(s1, -1),
  254.               decode8(s2, 0)));
  255.         ss->s1 = s1, ss->s2 = s2;
  256.     }    break;
  257.     case cEncode+cBits8+4:
  258.     {    register byte s1 = ss->s1, s2 = ss->s2, s3 = ss->s3;
  259.         loopn(4, (encode8(s0, -3), encode8(s1, -2),
  260.               encode8(s2, -1), encode8(s3, 0)));
  261.         ss->s1 = s1, ss->s2 = s2, ss->s3 = s3;
  262.     }    break;
  263.     case cDecode+cBits8+4:
  264.     {    register byte s1 = ss->s1, s2 = ss->s2, s3 = ss->s3;
  265.         loopn(4, (decode8(s0, -3), decode8(s1, -2),
  266.               decode8(s2, -1), decode8(s3, 0)));
  267.         ss->s1 = s1, ss->s2 = s2, ss->s3 = s3;
  268.     }    break;
  269.  
  270. #undef encode8
  271. #undef decode8
  272.  
  273.     }
  274. #undef loopn
  275. #undef dloop1
  276.     ss->row_left += count;        /* leftover bytes are possible */
  277.     if ( ss->row_left == 0 )
  278.     {    if ( end_mask != 0 )
  279.           *q = (*q & ~end_mask) | (save_last & end_mask);
  280.         if ( p < pr->limit && q < pw->limit )
  281.           goto row;
  282.     }
  283.     ss->s0 = s0;
  284.     pr->ptr = p;
  285.     pw->ptr = q;
  286.     return status;
  287. }
  288.  
  289. /* Stream templates */
  290. const stream_template s_PDiffE_template =
  291. {    &st_PDiff_state, s_PDiffE_init, s_PDiff_process, 1, 1, NULL,
  292.     s_PDiff_set_defaults, s_PDiff_reinit
  293. };
  294. const stream_template s_PDiffD_template =
  295. {    &st_PDiff_state, s_PDiffD_init, s_PDiff_process, 1, 1, NULL,
  296.     s_PDiff_set_defaults, s_PDiff_reinit
  297. };
  298.